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Abstract 

We derive a novel chiral power counting scheme for in-medium chiral perturbation theory with 
explicit nucleonic and pionic degrees of freedom coupled to external sources. It allows for a systematic 
expansion taking into account local as well as pion-mediated inter- nucleon interactions. Based on this 
power counting, one can identify classes of non-perturbative diagrams that require a resummation. 
Within this scheme, the pion self-energy in asymmetric nuclear matter is analyzed and calculated up- 
to-and-including next-to- leading order (NLO). It is shown that the corrections involving in- medium 
nucleon- nucleon interactions cancel between each other at NLO. As a result, there are no corrections 
up to this order to the linear density approximation for the in-medium pion self-energy. 



1 Introduction 



One of the long-standing issues in nuclear physics is the calculation of atomic nuclei and nuclear matter 
properties from microscopic inter-nucleon forces in a systematic and controlled way. This is a non- 
perturbative problem involving the strong interactions. In the last decades, Effective Field Theory 
(EFT) has proven to be an indispensable tool to accomplish such an ambitious goal. The pions play 
a unique role in the physics of the strong interactions. They should be included consistently with the 
spontaneous symmetry breaking of the SU(2)i ® SU(2)r chiral symmetry of the strong interactions in 
the limit of massless u and d quarks. We follow the techniques of Chiral Perturbation Theory (CHPT) 
[H El [3], with nucleons and pions as the pertinent degrees of freedom. In this approach, the pions are 
the Goldstone bosons associated with the spontaneous symmetry breaking of the chiral symmetry. They 
finally acquire a finite mass because of the small but non-vanishing masses of the lightest quarks u and 
d, which explicitly break chiral symmetry. For the lightest nuclear systems with two, three and four 
nucleons, CHPT in nuclear systems has been successfully applied [U EJ El EJ [9j [10] . For heavier nuclei 
one common procedure is to employ the chiral nucleon-nucleon potential derived in CHPT combined with 
standard many-body methods, sometimes supplied with renormalization group techniques [TT]. One of 
the most pressing issues of interest is the consistent inclusion of multi-nucleon interactions involving 
three or more nucleons in nuclear matter and nuclei, see e.g. |12|. [T3l [T4l [15]. 

In ref.[16] many-body field theory was derived from quantum field theory by considering nuclear 
matter as a finite density system of free nucleons at asymptotic times. Based on these results ref.|17] 



derived a chiral power counting in the nuclear mediumJ^J However, only nucleon interactions due to pion 
exchanges were considered. In this work we derive an extended (and newly organized) power counting 
that takes into account local multi-nucleon interactions simultaneously to the pion-nucleon interactions. 
Many present applications of CHPT to nuclei and nuclear matter [131 HH [20l [HI [T71 [22l [23J El], only 
consider meson-baryon chiral Lagrangians (see e.g. [10] for a summary), without constraints from free 
nucleon-nucleon scattering. In addition, as it is well known since the seminal papers of Weinberg [2j [3] , 
the nucleon propagators do often count as the inverse of a nucleon kinetic energy 0(p~ 2 ), so that they 
are much larger than assumed. This, of course, invalidates the straightforward application to nuclear 
physics of the pion-nucleon power counting valid in vacuum, as used e.g. in refs. [25j [191 [20J [21]. 

All these calculations share the assumption that spontaneous chiral symmetry breaking still holds 
for finite density nuclear systems. This assumption can be cross-checked by calculating the tempo- 
ral pion decay constant in the nuclear medium [17J. Namely, let us consider the axial- vector cur- 
rent A 1 ^ = g(x)7^75(r*/2)g(x), with q{x) a two-dimensional vector corresponding to the quarks fields 
and t % the Pauli matrices. Spontaneous chiral symmetry breaking results because the axial charge, 
Q l A = J cPxAq(x), does not annihilate the ground state, denoted by \Vt). As long as the matrix element 
(Q\Q l A \TT a (p)) = z(27r) 3 £(p)po ft is not zero, the vacuum is not left invariant by the action of the axial 
charge and spontaneous chiral symmetry breaking happens. In the previous equation |7r a (p)) denotes a 
pion state with Cartesian coordinate a, three-momentum p, energy pq and ft is the temporal weak pion 



decay coupling[*j Note that due to the presence of the nuclear medium one should distinguish between 
the spatial and temporal couplings of the pion to the axial- vector current. The calculations in ref.|17] 
indicate a linear decreasing of ft with density, ft = /^(l — (0.26 ± 0.04)/o/ / oo), where f w = 92.4 MeV is 
the weak pion decay constant in vacuum and po is the nuclear matter saturation density. This result 
clearly indicates that it makes sense to use chiral Lagrangians in the nuclear medium up to central nu- 



#1 This approach was later generalized to finite nuclei and e.g. applied to the calculation of the pion-nucleus optical 
potential |18| . 

#2 Mathematically one can obtain meaningful results from i8(p)po in the chiral limit by considering wave packets [17] . 
/ dp|pr 1 /(p 2 )K a (p)) with /(0)=constant |2B]. 
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clear densities. Nonetheless, a more thorough calculation of ft within our present approach, including 
short range nucleon-nucleon correlations, should be pursued in order to check whether the dependence in 
density of ft remains stable or is subject to significant corrections. On the other hand, the form of the 
chiral Lagrangians changes depending whether the quark condensate {£l\uu + dd\£l) is large or small. In 
the former case we have standard CHPT [1] [27] and in the latter the so called generalized CHPT would 
result [28J. In SU(2) CHPT it has been shown that the first case holds [29]. Refs. pH [30] obtain that 
the in-medium quark condensate decreases linearly with density as 1 — (0.35 ± 0.09) p/po [IZ]- Thus, the 
standard CHPT scenario holds up to nuclear matter saturation density. Ref. [31] calculated higher order 
corrections to this result and found the same linear trend for symmetric nuclear matter up to p ~ po- 
For higher densities, the linear decreasing is softened and frozen to a reduction of a 40% with respect 
to the vacuum value. For the pure neutron matter the higher order corrections calculated in [31] do not 
spoil the linear decrease of the quark condensate even for densities p > po. 

Our novel power counting is applied to the problem of calculating the pion self-energy in asymmetric 
nuclear matter at next-to- leading order (NLO). This problem is tightly connected with that of pionic 
atoms [32|, [33] due to the relation between the pion self-energy and the pion-nucleus optical potential. 
Despite being an old subject, a conclusive calculation of the pion self-energy in a systematic and controlled 
expansion is still lacking. For recent calculations see |TTl [3U [23], [35], [25j [36] . In particular, the problem 
of the missing S-wave repulsion, the renormalization of the isovector scattering length a~ in the medium 
[Ml [23] and the energy-dependence of the isovector amplitude [33] is not settled yet, despite the recent 
progresses [37J [331 E] ■ 

Relativistic field theories of nuclear phenomena featuring manifest Lorentz covariance are widely used 
to describe properties of nuclear matter and finite nuclei. Typically nucleons are described as point Dirac- 
particles moving in large isoscalar scalar and vector mean fields generated self-consistently [38] . The scalar 
mean field drives a strong attraction of order 300-400 MeV at nuclear density and almost and equally 
strong vector repulsion. This is a benchmark characteristic of the so called Quantum Hadrodynamics 
|39[ [40l |4T| 142] . As discussed in ref. [40], when the empirical low-energy nucleon-nucleon scattering is 
described in a Lorentz-covariant fashion, it contains strong scalar and four- vector amplitudes [33]. These 
important aspects are kept in the previous models which are widely applied for studying nuclear matter 
and nuclei. Applications are calculated with different degrees of refinement since the aio mean-field model 
of [38], including more vector and scalar fields and additional renormalizable scalar meson self-couplings 
|40] with adjustable parameters. In our approach the important nucleon-nucleon dynamics is generated 
by applying chiral EFT to systems with nucleons and pions, as commented above. No explicit mean fields 
are included, but this is not at odds with the mean-field models. Our approach is a dynamical one that 
should reproduce the physical effects of such mean fields in terms of the self-interactions of the explicit 
degrees of freedom included. 

After this introduction, we derive in section [2] a novel chiral power counting in the medium that 
takes into account multi-nucleon local interactions, pion exchanges and the enhancement of nucleon 
propagators. In sections [3] and [3] we calculate the different contributions to the pion self-energy in 
asymmetric nuclear matter up-to-and-including 0(p 5 ) at NLO. Interestingly, we show in section [4] that 
the different NLO contributions from nucleon-nucleon scattering cancel between each other. Conclusions 
are given in section [5j 

2 Chiral Power Counting 

In ref.pjBj the effective chiral pion Lagrangian was determined in the nuclear medium in the presence of 
external sources. For that the Fermi seas of protons and neutrons were integrated out making use of 
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Figure 1: Representation of an IGV. See the text for further details. 

functional techniques. A similar approach was followed in ref.[H] but for the case of only one nucleon. In 
this way it is manifestly shown that pion or nucleon field redefinitions do not affect physical observables 
also in nuclear matter because they appear as integration variables in a functional. Nonetheless, in 
ref.[16] only the meson-baryon chiral Lagrangian is employed. More precisely, if we write a general chiral 
Lagrangian in terms of an increasing number of nucleon fields ip, 

A?ff = ArTr + Cfo/f + + • • • (2-1) 

only the contributions from C n7T and were retained in ref . [16j . Based on these results, the authors 
of ref.[T7j derived a chiral power counting in the nuclear medium. 

Ref. [116] establishes the concept of an "in- medium generalized vertex" (IGV). Such type of vertices 
result because one can connect several bilinear vacuum vertices through the exchange of baryon propa- 
gators with the flow through the loop of one unit of baryon number, contributed by the nucleon Fermi 
seas. This is schematically shown in fig[TJ where the thick arc segment indicates an insertion of a Fermi 
sea. At least one is needed because otherwise we would have a vacuum closed nucleon loop that in a low 
energy effective field theory is buried in the higher order chiral counterterms. On the other hand, a filled 
large circle in figH] indicates a bilinear nucleon vertex from £ n N, while the dots refer to the insertion 



of any number of them|*f| It was also stressed in ref. 1 17] that within a nuclear environment a nucleon 
propagator could have a "standard" or "non-standard" chiral counting. To see this note that a soft 
momentum Q ~ p, related to pions or external sources attached to the bilinear vertices in figdj can be 
associated to any of the vertices. Denoting by k the on-shell four-momenta associated with one Fermi 
sea insertion in the IGV, the four-momentum running through the j th nucleon propagator can be written 
as pj = k + Qj . In this way, 



{k + Qj) 2 -m 2 + ie Q 2 + 2Q° j E(k) - 2Qjk + ie ' 

where E(k) = k 2 /2m, with m the physical nucleon mass (not the bare one), and Q® is the temporal 
component of Qj. We have just shown in the previous equation the free part of an in-medium nucleon 
propagator because this is enough for our present discussion. Two different situations occur depending 
on the value of Q®. If Q® = 0(m n ) = 0(p) one has the standard counting so that the chiral expansion 
of the propagator in eq. H2.2p is 



2Q v j m + ie V 2Q {J j m 

Thus, the baryon propagator counts as a quantity of 0(p~ l ). But it could also occur that Q® is of the 
order of a kinetic nucleon energy in the nuclear medium or that it even vanishes. The dominant term in 



#3 /•(*) 



£jriv corresponds to the CHPT pion-nucleon Lagrangian of chiral order i. 
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eq. (jlT2]) is then 



y, + Qj +m 



(2.4) 



— i 



Q2 + 2Qj • k - is 



and the nucleon propagator should be counted as 0(p~ 2 ), instead of the previous 0(p^ 1 ). This is 
referred to as the "non-standard" case in ref . [TTj . We should stress that this situation also occurs already 
in the vacuum when considering the two-nucleon reducible diagrams in nucleon-nucleon scattering. This 
is indeed the reason advocated in ref. [2] for solving a Lippmann-Schwinger equation with the nucleon- 
nucleon potential given by the two-nucleon irreducible diagrams. In the present investigation, we extend 
the results of refs. |16|. \T7\ in a twofold way. i) We are able to consider chiral Lagrangians with an 
arbitrary number of baryon fields (bilinear, quartic, etc). First only bilinear vertices like in refs. |16^ 117] 
are considered, but now the additional exchanges of heavy meson fields of any type are allowed. The 
latter should be considered as merely auxiliary fields that allow one to find a tractable representation 
of the multi-nucleon interactions that result when the masses of the heavy mesons tend to infinitvl^l 
These heavy meson fields are denoted in the following by H, see figGJ and a heavy meson propagator is 
counted as O(p ) due to their large masses, ii) We take the non-standard counting from the start and 
any nucleon propagator is considered as 0(p~ 2 ). In this way, no diagram whose chiral order is actually 
lower than expected if the nucleon propagators were counted assuming the standard rules is lost. This is 
a novelty in the literature. 

In the following ~ kp ~ 0{p) are taken of the same chiral order, and are considered much smaller 
than a hadronic scale A x of several hundreds of MeV that results by integrating out all other particle 
types, including nucleons with larger three-momentum, heavy mesons and nucleon isobars [3]. The chiral 
order of a given diagram is represented by v and is given by 



From left to right, Lh is the number of loops due to the internal heavy mesonic lines, L n that of pionic 
loops and I n is the number of internal pionic lines. Each loop introduces a factor of four in the power 
counting, because of the integration over the free four-momentum, and a pion propagator reduces the 
order by two units. The quantity di is the chiral order of the i th bilinear vertex in the baryonic fields 
and V is the total number of such vertices. V p is the number of IGVs and mi is the number of nucleon 
propagators in i th IGV minus one, and every of them reduces by two units the chiral counting as discussed 
previously. The definition of rrii contains the removal of one baryon propagator because there is always 
at least one Fermi sea insertion for each IGV that increases the chiral counting in three units because 
the associated integration over a Fermi sea, f d?k9{^i — |k|), with £j the corresponding Fermi momentum. 
This fact gives rise to the last sum in eq . (j2.5j) . Other symbols that appear are £i, that is the chiral 
order of a vertex without baryons (only pions and external sources), and which is the total number 
of the latter ones. We have not included in eq. (|2.5p any contribution from tt-H vertices without baryons 
because in the limit when the mass of the H fields is taken to infinity the H propagators are contracted 
to a point and the pions will be always attached to baryons. 

Let us note that associated with the bilinear vertices in an IGV one has four-momentum conservation 
Dirac delta functions that can be used to fix the momentum of each of the baryonic lines joining them, 
except one for the running three-momentum due to the Fermi sea insertion. Let us now introduce another 
symbol, V$. Here, we take as a whole any set of IGVs that are joined through heavy mesonic lines H, 

* 4 Such methods are also used in the so-called nuclear lattice simulations, see e.g. [45] . 




V = 



(2.5) 



i=l i=l i=l i=l 
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Figure 3: Representation of a possible arrangement of IGVs separated in two clusters. In this figure V p = 5, 
V$ = 2, ijr = 3, Ih = 3 and E = 1. The pions are indicated by the dashed lines and the external source by a wavy 
line. Eqs. (|2.7p and (|2.6[) imply that L n = 2 and Lh = as it should. 



whose total number is Ih- The number of these clusters of in-medium generalized vertices is denoted by 
V$. In this way, we can write 

L H = I H ~Y J {V P ,i - 1) = Ih ~ V p + V* , (2.6) 

i=l 

where V P j is the number of IGVs within the i th set of generalized vertices connected by heavy mesonic 
lines. Additionally, they could be connected between them or with other IGVs belonging to other clusters 
by pionic lines. Since there is a total four-momentum conservation delta function associated to every of 
these clusters it follows that 

L v = J ff - V v - + 1 . (2.7) 

These relations are illustrated in figj3] where a possible arrangement of IGVs is shown. On the other 
hand, 

V 14 

2I H + 2I^ + E = J2 v i + J2 ni ■ ( 2 - 8 ) 

i=l i=l 

Here, V{ is the number of mesonic lines attached to the i th bilinear vertex, n« is the number of pions in 
the i th mesonic vertex and E is the number of external pions. Taking into account eqs. (|2.6p . (|2.7p and 
eq-CEH]), ea. ([23]l reads , 

V„ V 

v = 2I H ~ E + 4 - 4K + Y^li + m) + Y^Ai + Vi) -2m-V p , (2.9) 

i=i i=i 

with m = Y2i=i m i ■ We now employ in eq. (USD that V p + m = V , and 21 H = Y%=1 ^ ' where 

LUi is the 

number of heavy meson internal lines for the i th bilinear vertex. Then, we arrive at our final equation 

V-n V m V p 

„ = 4 - E + Y J (n l + £i - 4) + Y(di + - 1) + Y J (vi ~ 1) + Y Vi ■ ( 2 - 10 ) 

i=l i=l i=l i=l 
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Note that v given in eq. (|2,10p is bounded from below because n.j + 1\ — 4 > 0, as li > 2 and n% > 2, 
except for a finite number of terms that could contain only one pion line but always having external 
sources attached to them. Similarly di + uj{ — 1 > 0. For the pion-nucleon Lagrangians this is always 
true as di > 1. For those bilinear vertices mediated by heavy lines di > but then Ui > 1. For the 
term before the last one in eq. (j2.1Qp vi — 1 > 0, except for the higher-order nucleon-mass renormalization 
counter terms or the finite number of terms which would not have pionic lines but only external sources 
from £ n N- The former terms have di > 2 and then (di + w$ — 1) + (u$ — 1) > 0. For di = 2 the chiral 
order does not increase but these terms can be absorbed in the physical nucleon mass. For the last term 
in eq. (l2,10p > and then it is positive. It is worth stressing that adding a new IGV to a connected 
diagram increases the counting at least by one unit because then V{ > 1. Using again that V p + m = V 
eq. (|2.10p can be rewritten as 

V n V 

v=4-E+ Y^irii +£i-4) + Y^{di + u>i + Vi - 2) + V p . (2.11) 

i=l »=1 

The number z/ given in eq. (|2.10p represents a lower bound for the actual chiral power of a diagram, 
f/,, so that /x > v. The real chiral order of a diagram might be different from v because the nucleon 
propagators are counted always as 0(p~ 2 ) to obtain eq. (|2.10p . while for some diagrams there could be 
propagators that follow the standard counting. Eq. (|2,10|) implies the following conditions for augmenting 
the number of lines in a diagram without increasing the chiral power by: 

1. adding pionic lines attached to mesonic vertices, t% = rii = 2. 

2. adding pionic lines attached to meson-baryon vertices, di = V{ = 1. 

3. adding heavy mesonic lines attached to bilinear vertices, = 0, = 1. 



There is no way to decrease the order We apply eq. (l2,10p by increasing step by step V p up to the 
order considered. For each V p we look for those diagrams that do not increase the order according to 
the previous list. Some of these diagrams are indeed of higher order and one can refrain from calculating 
them by establishing which of the nucleon propagators scale as 0(p^ 1 ). 



3 Meson-baryon contributions to the pion self-energy 

Let us apply the chiral counting given in eq. (|2.1ip to calculate the pion self-energy in the nuclear medium 
up to NLO or 0(p 5 ), with the different contributions shown in figHl The nucleon propagator, Go(k)i 3 , 
contains both the free and the in-medium piece 



In this equation the subscript i% refers to the type of nucleon, so that, = p corresponds to the proton 
and 13 = n to the neutron. Our convention for the pion self-energy, E, is such that the dressed pion 
propagator reads 

qz _ m z _|_ 



#5 Only by adding vertices with li = 2 and rn < 2 or di — 1 and Vi — 0. However, its number is bounded from above by 
the necessarily finite number of external sources. 
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o( P 4 ) 

Leading Order 



V p = l 
0(p 5 ) 

Next-to-Leading Order 



V p = 2 

Next-to-Leading Order 



Q 

f 3 










X ~>~<r >-< x: 



Figure 4: Contributions to the in- medium pion self-energy up to NLO or 0(p 5 ). The pions are indicated by 
the dashed lines and the squares correspond to NLO pion-nucleon vertices. The wiggly lines correspond to the 
nucleon-nucleon interaction kernel, see figE^,, that is iterated as meant by the ellipsis. 



We employ the pion-nucleon Heavy Baryon CHPT (HBCHPT) Lagrangians at 0(p) and 0(p 2 ) 
can be found e.g. in refs.|47t 148] . For completeness we reproduce them, 



that 



C 



(i) 

ttN 



N(iD 



9A 



a -u)N 



C 



(2) 
rJV 



N 



( 2i? 5 ' 3 + ij m {° ' 3 > u °} + 2cim *( u + u ^ + ( C2 - w) n ° + csU » uti ) + • • • ( 3 - 3 ) 



where the ellipses represent terms that are not needed here. In this equation, N is the two component field 
of the nucleons, qa is the axial pion-nucleon coupling and = + the covariant chiral derivative, 
being = [u',dn,u]. The pion fields tv(x) enter in the matrix u = exp(ir • 7f/2/), in terms of which 
Uf X = i {u^c^ii} and U = u 2 , with / the weak pion decay constant in the SU(2) chiral limit. The c, are 
chiral low energy constants whose values are fitted from phenomenology [48J. 

The leading contribution to the pion self-energy corresponds to the diagrams 1-3 on figiH The diagram 
1 results by closing the Weinberg- Tomozawa pion-nucleon interaction (WT), obtained by expanding £^ 
in eq. (|3.3p up to two pion fields, 



V' 



-£ij3(Pp — Pn) 



(3.4) 



8 



where the proton(neutron) density is given by p p ( n ) = £p( n )/37r 2 . Si is then a S-wave isovector self-energy. 
The sum of the diagrams 2 and 3 of figH] is 

ig 2 q 2 g 2 (q 2 ) 2 

S 2 = -^0 £ ij3(Pp ~ Pn) ~ tp-^r^jiPp + Pn) ■ (3-5) 

They involve the one-pion vertex from C^lL, which is proportional to gA- This is a P-wave self-energy 
where the first term is isovector and the second is isoscalar. 

Now, we move to the NLO contributions. The vertices from C ^ in eq. (|3.3p are indicated by squares 
in the figlH It should be understood that the pion lines can leave or enter the diagrams 4 and 5 of figlH 
The sum of these two diagrams gives the result, 

X3 = ^(Pp + Pn)S ij • (3.6) 

(2) 

This is a P-wave isoscalar contribution that stems from the term linear in g^ in C^L, that is a recoil 
correction to . The diagram 6 of figlH is given by 

2 



P 



' 2 ' ),j f 2ci ml - q 2 (c 2 + c 3 - J^) + c 3 q 2 ) (p p + p n ) , (3.7) 



in terms of the q couplings of C^Jf, eq. (|3.3|) . This is an isoscalar contribution where the term — 25ijC^q 2 (p p + 
Pn)/ f 2 is P-wave and the rest is S-wave. 

Next, let us consider the contributions to the pion self-energy due to the nucleon self-energy from a 
one-pion loop as depicted in the diagrams 7-9 of figSJ with vertices calculated from The diagrams 

originate by dressing the in-medium nucleon propagator of the diagrams 1-3 by the one-pion loop nucleon 
self-energy, 

1+T3^ + 1_T3 S ^ (38) 

with and the proton and nucleon self-energies due to the in- medium pion- nucleon loop. Their 
values at k° = are subtracted because we employ the physical nucleon mass. The contribution from 
the diagram 7 of figH] can be written as 

S 5 = |U' 3 / (04 - G (k) 2 n K) ^ , (3.9) 

where e zk ° v , rj —* + , is the convergence factor associated to any closed loop made up by a single nucleon 
line [46]. Taking into account that 



G (k) 2 3 = -dGo{k) i3 /dk° , (3.10) 
we then integrate by parts, which is possible thanks to the convergence factor. It results, 

Ss = J 2 ^ J [ Go{k) ^ ~ Go{k)n m£) e " ■ (3 - n) 



Since r\ — > the additional term obtained by taking the derivative of e lfc0r? with respect to k° in eq. (|3.9p 
does not contribute and is not shown. £5 is an isovector S-wave pion self-energy contribution. For the 
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diagrams 8 and 9 of the same figure one has analogously 



o fcl i 3 



/ 



•2 



d 4 k 

d 4 k , 



G (k) 



n dk° 



ik°rj 



G (fc) p S^ + Go(fe) n E^)e 



(3.12) 



Eg is a P-wave self-energy contribution but while the first line is of isovector character, the one in the 
second line is isoscalar. This last term is indeed a NNLO or C(p 6 ) contribution because the pion-loop 
nucleon self-energy is 0(p 3 ) and we neglect it. The free pion-loop nucleon self-energy is calculated in 
HBCHPT |48| . employing dimensional regularization in the MS — 1 scheme. Its derivative is 



dT, l{n)J 
dk° 



32vr 2 / 2 



+ k, 



3k°Vb \ Hog 



k° + iVb 
-k° + iVb 



(3.13) 



with b = m\ — k\ — ir\. Hence, because dTi^^ j./dk° = 0{p 2 ) when inserted in £5 and T,q it gives rise to 
an 0(p e ) contribution that we neglect in the present work. The in-medium contribution to the pion-loop 
nucleon self-energy involves the finite integral 



2tt 



dH es(k°-i°)e(c 



'3 



(27^ 



mi + ie 



dH P0(£. 



(2tt) 3 b + P 



(3.14) 



This integral only depends on k° through the variable b = m\ — k\. Since dl m /dk° = —2k°dl m /db = 
0(p 2 ), because k° = 0(p), the in- medium part of the pion loop contribution to the nucleon self-energy 
gives rise to 0(p e ) terms for the pion self-energy. As a result, £5 and £6 are at least 0(p G ). Notice 
that from eq. (|2.10p these contributions were firstly booked as (D(p 5 ) because it was considered that 
9£ 7r /5fc° = 0(p) as £ 7r = 0(p 3 ) and k° = 0{p 2 ). But since for these diagrams V p = 1, with only one 
closing nucleon line, k° = 0(p). 



4 In-medium nucleon-nucleon scattering contributions 

We now consider those NLO contributions to the pion self-energy in the nuclear medium that involve the 
nucleon-nucleon interactions. They are depicted in the diagrams of the last two rows of figSl where the 
ellipsis indicate the iteration of the two- nucleon reducible loops. For the diagrams b) and d) of fig H] the 
pion lines can leave or enter the diagrams. It is remarkable that these NLO contributions cancel between 
each other. On the other hand, since V p = 2 in these contributions one needs only the nucleon-nucleon 
scattering amplitude at O(p ) to match with our required precision at NLO. This amplitude is obtained 
by iterating in an infinite ladder of two nucleon reducible loops, with full in-medium nucleon propagators, 
the tree level amplitudes obtained from the O(p ) Lagrangian with four nucleons [3] and from the one- 
pion exchange with the lowest order pion-nucleon coupling. This procedure would correspond in vacuum 
to the leading nucleon-nucleon scattering amplitude according to refs.[2j[3]. 
The O(p ) lowest order four nucleon Lagrangian [3] is 

£ ( nn = -\Cs(NN)(NN) - ±C T (NaN)(NaN) . (4.1) 

Of course, this Lagrangian only contributes to the S-wave nucleon-nucleon scattering. The scattering 
amplitude for the process N slj i l (pi)N S2 ^ 2 (p2) — ► A r S3i j 3 (p3)A r S4j j 4 (p4), with s m a spin label and i m an 
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a) b) 



Figure 5: a) The wiggly line represents the sum of the leading one-pion exchange and four-nucleon contact inter- 
actions, T^ N + T}f N , eqs. (|4.2p and (|4.3|) . b) In-medium nucleon self-energy due to the nucleon-nucleon interactions 
with the Fermi seas. 



isospin one, that follows from the previous Lagrangian is 

^NN — {^S3Si^S4S2 ^1311^1^2 ^S3S2^S4Si ^1342^*4*1) 

C7 1 ' &S4S2 ^1311^1412 ^S3S2 ' ^1312*^1411) - (^"^) 

Regarding the one-pion exchange at leading order, with vertices from one has the expression 



-,2 



in _ 9 a (%i • ^4i 2 )(^ • q)s 3 si(o ! ' <l)*i«a (ruh ■ T i3 i 2 ){B ■ q!) S4S1 (a • q') S3 s 2 



T NN 4f 



q 2 + m 2 - ie q' 2 + m 2 - ie 

PC I rplTT 

■ NN 1 NN 



(4.3) 



with q = P3 — Pi and q' = P4 — Pi- In the following, the sum T^ N + T^ N is represented diagrammatically 
by the exchange of a wiggly line, as in figl5^. 

The diagrams a) and c) of figH] involve the WT vertex while b) and d) contain the pole terms of 
pion-nucleon scattering. At leading order in the chiral counting the sum of the latter two has the same 
structure as the WT term, with the resulting vertex given by 

^-(l-gl^)e ijk r k . (4.4) 



2/ 2 V ^ 

We can then discuss simultaneously all the diagrams in the last two rows of figlU employing the latter 
vertex. The sum of the diagrams a) and b) of fig|4] can be written in terms of the nucleon self-energy in 
the nuclear medium due to the nucleon-nucleon scattering. It reads 

£7 = ^ (l - 9 2 A^) e ijk J l^e**" Tr jr* (^G (k l)p + ^<*,(fci)„) V NN (^G (k 2 ) P 
+^P L G (k 2 ) n )\ . (4.5) 



Here, 



^NN = ^— -^p,NN H 7T~ ~^n,NN , (4.6) 



with Sp( n ) iA rAf the proton (neutron) self-energy in the nuclear medium due to the nucleon-nucleon inter- 
actions, fig Eh, 

/d^k o 
jY\4 G o( k 2)a2TNN(kicrii 3 ,k 2 a2a2\k 1 a 1 i 3 ,k2cr2a2)e tk ° v . (4.7) 

In this expression Tnn is the nucleon-nucleon scattering amplitude between the indicated initial and 
final off-shell states. These are characterized by three labels. The first label corresponds to the four- 
momentum, the second to the spin and the third to the isospin. Note that in the equation there is a sum 
over all the quantum numbers of the second nucleon. 
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a) b) 

Figure 6: Contribution to the pion self-energy with a two-nucleon reducible loop. The pion scatters inside/outside 
the loop for the diagram a) /b) . 



Proceeding similarly as before for £5, including the integration by parts, we can write 



Inserting in the previous equation the explicit expression for E^aw of eq. (|4.7p one has 



a 



d A h d 4 k 2 ik o v ik o v 
(2^) 4 (2^) 4 



G (ki) p G (k2)p-g^T NN (kiaip, k 2 a 2 p\k l a 1 p, k 2 a 2 p) 
d 

G (k 1 ) n G (k 2 ) n -^T NN (k 1 a 1 n, k 2 a 2 n\k 1 ain, k 2 a 2 n) 



This result is obtained from eq. (|4.8p by noting that 



d 4 h d A k 2 _ e ik°r, e ik° V 



(2vr) 4 (2vt) 4 



d 



Go(k 1 ) p G (k2)n-Q^T NN (k 1 a 1 p, k 2 a 2 n\kxaxp, k 2 a 2 n) 
d 

- G (k 1 ) n G {k 2 )p-^T NN (k 1 a 1 n, k 2 a 2 p\k 1 a 1 n, k 2 a 2 p) 

This follows for two reasons. First, let us notice that because of Fermi-Dirac statistics 
T NN (k 1 aip,k 2 a 2 n\kia 1 p,k 2 cr 2 n) = T NN (k 2 a 2 n,kiaip\k 2 a 2 n,k 1 aip) . 



(4.8) 



(4.9) 



(4.10) 



(4.11) 



Second, at LO the amplitude T/vat, as commented above, is given by the iteration of the wiggly line in 
figEJ The latter does neither depend on k^ nor on fcf>, see eqs. (|4.2p and (|4.3p . Since the two-nucleon 
reducible unitarity loop, the one in ngEfr, depends on k® and k® through their sum, k^ + k®, as can be 
seen straightforwardly [49J, it results that T/vat at LO only depends on them in this way. It follows then 
that BTnn / dk® = OTnn /dk®. Taking these two facts into account, as ki and <7j are dummy variables, 
eq. fl4.10p results. 

Let us now consider the diagrams c) and d) of figH] whose contribution is denoted by Eg. These 
diagrams consist of the pion-nucleon scattering in a two-nucleon reducible loop which is corrected by 
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initial and final state interactions. The iterations are indicated by the ellipsis on both sides of the 
diagrams. In order to see that these diagrams cancel with eq. (|4.9p let us take first the diagram of figl6^ 
with a twice iterated wiggly line vertex. It is given by 



d 4 h d A k 2 

(2vr) 4 (2vr) 4 



G (k 1 ) a G (k 2 ) l 3e ik ° T 'e ik ° r ' 



a', 13' a[,a' 2 



dG (ki - q) a > 



T 3\a'a' V a'(3';ap{-q)Go{k 2 + q)p 



(4.12) 



where V a j3- n $ corresponds to the wiggly line with the indices a and 7 belonging to the out-/in-going 
first particle, in that order, and similarly (3 and 5 for the second one. To shorten the notation, we have 
only indicated the isospin indices in V in the previous equation, although V depends also on spin. A 
symmetry factor 1/2 is also included because V a p^g contains both the direct and exchange terms, as 
explicitly shown in eqs, (|4.2p and (14, 3D . The derivative with respect to k® arises in eq, (!4,12j ) because 
the nucleon propagator to which the two pions are attached appears squared and we have made use 
of eq. (13.10p . This is so because for the ir^, i and j can be either 1 or 2, so that the only surviving 
contribution is k = 3. For the 7r°, i = j = 3 and then there is no contribution. Thus, because one 
has either or r 3 , which is a diagonal matrix, the nucleon propagator before and after the two-pion 
vertex is the same. The q-loop integral in eq. (|4.12p is typically divergent. Nevertheless, the parametric 
derivative with respect to k® can be extracted out of the integral as soon as it is regularized. Of course, 
the same regularization method as that used to calculate T/vjv should be employed. Once the derivative 
is taken out of the integral, the quantity between the squared brackets in eq. (|4.12p corresponds to the 
twice iterated wiggly line contribution to TNN(k%a%a, k 2 cr2l3\kiaia, k2(T%0). In addition, the isovector 
nature of the modified WT vertex of eq. (|4.4p implies that only the difference between the proton-proton 
and neutron- neutron contributions arises. This can be worked out straightforwardly from the isospin 
structure of the local four-nucleon vertex and that of the one-pion exchange [49j. As a result we can 
write 



d A ki d 4 k 2 
(2vr) 4 (2tt) 4 







d 

G (k 1 ) p Go(k 2 )p-^T^ N (k 1 a 1 p, k 2 a 2 p\k x o x V^ k 2 a 2 p) 



G (ki) n G (k2)n-7^-pTfc N (ki(Tin, /c 2 o- 2 n|£;i(7in, k 2 a 2 n) 

where the superscript L on top of Tjvat indicates that this scattering amplitude is calculated at the one- 
loop level. This result then cancels exactly with that of figEb, corresponding to the twice iterated wiggly 
line exchange contribution to Tnn in eq. (|4.9|) . This cancellation is explicit by reducing eq. (|4.9p to the 
one- loop case for calculating Tjvtv Notice as well that the contribution to TV at given by the exchange of 
only one wiggly line vanishes when inserted in eq. (|4.9p because it is independent of see eqs. (|4.2p and 

(HP- 

This process of mutual cancellation between S7 and Ss can be generalized to any number of two- 
nucleon reducible loops in figsIUa), b) andHfc) and d), respectively. An n + 1 iterated wiggly line exchange 
in these figures implies n two-nucleon reducible loops. The two pions can be attached for Ss to any of 
them, while for £7 the derivative with respect to can also act on any of the loops. The iterative loops 
are the same for both cases but a relative minus sign results from the loop on which the two pions are 
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Figure 7: In this figure the cross indicates the action of the derivative with respect to k° in eq. (|4.9p . When the 
derivative is performed over a baryon propagator the latter becomes squared, according to eq. (|3.10|) . In this way, 
the first diagram on the second row of the figure is the same as the one to the left of the first row but with opposite 
sign and they cancel each other. The same applies to the second diagrams on the first and second rows. 

attached with respect to the one on which the derivative is acting, as just discussed. This is exemplified 
in figJT] for the case with two two-nucleon reducible loops. Hence, 

S 7 + S 8 = . (4.14) 

The basic simple reason for such cancellation is that while for £7 the presence of a nucleon propagator 
squared gives rise to {-ifd/dk^, for S 8 it yields -d/dk%, cf. eqs.gji]) and IfOSjl . respectively. The 
extra (—1) for E7 results because it involves an integration by parts, as discussed above. 

5 Conclusions and outlook 

We have developed a promising scheme for an EFT in the nuclear medium that combines both short- 
range and pion-mediated inter- nucleon interactions. It is based on the development of a new chiral power 
counting which is bounded from below and at a given order it requires to calculate a finite number 
of contributions. The latter could eventually involve infinite strings of two-nucleon reducible diagrams 
with the leading (D(p ) two-nucleon CHPT amplitudes. As a result, our power counting accounts for 
non-perturbative effects to be resummed which, e.g., give rise to the generation of the deuteron in 
vacuum nucleon-nucleon scattering. The power counting from the onset takes into account the presence 
of enhanced nucleon propagators and it can also be applied to multi-nucleon forces. 

We have then calculated the leading corrections to the lowest order result for the pion self-energy in 
asymmetric nuclear matter, with all the contributions up-to-and-including 0(p 5 ) evaluated. As a novelty, 
it is shown that the leading corrections to the linear density approximation vanish. In particular, it is 
derived that the leading corrections from nucleon-nucleon scattering mutually cancel. This suppression 
is interesting since it allows to understand from first principles the phenomenological success of fitting 
data on pionic atoms with only meson-baryon interactions [3T^ I33j . An 0(p G ) calculation of the pion self- 
energy is a very interesting task as it provides the first corrections to the linear density approximation, 
e.g. the well-known Ericson-Ericson-Pauli rescattering effect [32] . More calculations for other physical 
processes and higher orders are clearly needed to assess the realm of applicability of the present approach. 
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